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ABSTRACT 

Using a simple method to interpret the luminosity function of Lya emitters, we ex- 
plore properties of Lya emitters from 5.7 ^ z ^ 8.8 with various assumptions about 
metallicity and stellar mass spectra. We constrain a mass-to- "observed light" ratio, 
Mh/Lband, where Mh refers to the total mass of the host halo, and Lband refers to the 
observed luminosity of the source. For narrow-band surveys, Lband is simply related to 
the intrinsic Lya luminosity with a survival fraction of Lya photons, aesc- The mass- 
to- "bolometric light", Mh/LboU can also be deduced, once the metallicity and stellar 
mass spectrum are given. The inferred Mh/Lboi is more sensitive to metallicity than to 
the mass spectrum. We find the following constraints on a mass-to-light ratio of Lya 
emitters from 5.7 ^ ^ < 7: {Mh/ Lboi){aesce^^''y^ 21 - 38, 14 - 26, and 9-17 for 
Z = 0, 1/50, and 1 Zq, respectively, where e is the "duty cycle" of Lya emitters, and 
7 ~ 2 is a local slope of the cumulative luminosity function, A''(> L) oc L~'^, to which 
the current data are sensitive. Only weak lower limits are obtained for z = 8.8. There- 
fore, Lya emitters are consistent with either starburst galaxies [Mh/Lboi ^ 0.1 — 1) 
with a smaller Lya survival fraction, aescC^^''' ~ 0.01 — 0.05, or normal populations 
[Mh/Lboi ~ 10) if a good fraction of Lya photons survived, aescC^^^ ~ 0.5 — 1. We 
find no evidence for the end of reionization in the luminosity functions of Lya emitters 
discovered in the current Lya surveys, including recent discovery of one Lya emitter 
at 2 = 7. The data are consistent with no evolution of intrinsic properties of Lya 
emitters or neutral fraction in the intergalactic medium up to z = 7. No detection of 
sources aX z — 8.8 does not yield a significant constraint yet. We also show that the 
lack of detection at 2; = 8.8 does not rule out the high-2 galaxies being the origin of 
the excess near infrared background. 

Key words: cosmology: early Universe, observations, theory - infrared: galaxies - 
galaxies: high redshift 



1 INTRODUCTION 

What was the universe like at high redshifts? We are cur- 
rently entering a time where we can begin to observe this 
early time in the universe's life by looking for galaxies at 
redshifts above 6. 

There are several indications that stars, and hence 
galaxies, existed at this early time. We know that the 
universe was reionized early from observations such as 
the polarized light of the cosmic mi c rowav e background 
fealdarr iagal Il997 ': 'Kaplingha t et all l2003l: iKogut et all 



12003. : Page et all |2007: Spergel et al. 2007), the Gunn 



Peterson test towards q u asars ( e.g. Gunn &: Petersonl 
Fan et al.l bool I2001I |2002|. [20041 : iBecker et al.l 



1965 



2001 



Oh fc Furlanettdl2005l : lGotoll2006i rand a gamma-ray burst 
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(|Totani et al."2006') , and the tem perature of the intergalac- 
tic medium (iHui fc HaimanlboOSi ). 

In order to produce large scale reionization, an efficient 
and plentiful source of ultraviolet photons was needed. The 
first few generations of stars are very likely candidates, pro- 
ducing ultraviolet photons efficiently. The Lya forest shows 
that the universe was polluted with metals as ea rly as z ~ 6 
(Songaila 2006; Pottini ct al. 2003; Ryan- Webe r et al.ll2006l : 
ISimcoel200 6). indicating even earlier star formation and thus 
providing further support of early stars. In addition, a por- 
tion of the near infrared background could be the redshifted 
light from t he first sta rs, and may provide i nformation 
about them jSantos et al. 2002; Salvat erra fc Fe rrara 200j; 



Gooray fc Yoshida 



20071 ) 



2004 



Mii fc Totani 2005; Madau fc Silk 



20051: [Fernandez fc Komatsul l20od iKashlinskv et al.l i2005r 



It is therefore very likely that there is significant star 
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formation above 2 > 6. With the introduction of new, 
more powerful telescopes and deep field searches, an in- 
teresting question arises: do these first stars form galaxies 
that are bright enough to be seen today? Several deep field 
Lya s earches have been performed with the ISAAC on the 
VLT (|Willis fc CourbinlbOOSl i lWiUis et al.|[2006l: ICubv et al. ' 



2007 ). the Mayall Telescope at Ki tt Peak (iRhoads et al 
2004 ) , an d the Subaru telescope (iTaniguch i et al. _ 20051: 



Ive et all I2OO6I : iKashikawa et al] [2OO6 : Shimasaku et all 
20061 ). probing the universe at redshifts of six and higher. 

In this paper, we present a simple method to calculate 
the luminosity function of high-z galaxies, and compare this 
with the results of Lya searches to constrain properties of 
Lya emitters. In § [2l we explain our method to calculate 
the luminosity function of high-z galaxies with a single free 
parameter, a mass-to-observed light ratio, Mh/ Li,and - In §|3l 
we review the current generation high redshift galaxy sur- 
veys and their data. In § U we model stellar populations to 
obtain actual physical quantities about the galaxies that we 
observe. We compare our results to previous work in § [5] 
and conclude in § [S] We use the latest WMAP 3 parame- 
ters of = 0.0422, = 0.241, Oa = 0-759 ag = 0.761, 
h = 0.732, and n, = 0.958 jSpergel et al.ll2007h . 



2 A SIMPLE MODEL OF GALAXY COUNTS 

2.1 Justification for a simplified approach 

The simplest way to predict the cumulative luminosity func- 
tion of galaxies is to count the number of haloes available in 
the universe above a certain mass. 



N(> L) = V{z) 



dn 



dMh, 



(1) 



where V{z) is the survey volume, dn/dMh the comoving 
number density of haloes per unit mass range, and Mh the 
total mass of a halo. 

The cumulative number density of haloes, 

dMh dn/dMh, is shown in the bottom panel of 
Figure [T] In order to calculate the cumulative luminosity 
function, one may simply "stretch" the horizontal axis of 
this figure by dividing Mh by a suitable factor that converts 
the mass to luminosity: a mass-to-light ratio, Mh/L. 

This model is admittedly oversimplified, and is indeed 
simpler than what is already available in the literature. For 
example, one can stretch not only the horizontal axis (i.e., 
mass) , but also the vertical axis of the cumulative mass func- 
tion as J dMh dn/dMh e f dMh dn/dMh, wh ere e is often 
called a "duty cycle" (e.g.. iHaiman et al.ll2000l ). 

The vertical stretch would be required when the aver- 
age lifetime of Lya emitting galaxies, Tg, is shorter than 
the age of the universe, in which case the number count 
should be given by the time derivative of the mass func- 
tion, / dMh dn/dMh J dMh Tg{Mh) d^n/{dMhdt). The 
vertical stretch parameter is thus given approximately by 

6 ~ Tg / i univ . 

Since the statistical power of the current data is not yet 
strong enough to constrain both the horizontal and verti- 
cal stretch paramete rs simultaneously, these parameters ar e 
strongly degenerate (jPiikstra et al.ll2006l : [Stark et al.ll2007h . 
They are completely degenerate when A^(> L) follows a sin- 



gle power law, A'^(> L) oc L~'^ . The degeneracy line is given 
by [Mh/ L)e~^^^ = constant; thus, the inferred Mh/L and 
e are positively correlated: the smaller the e is, the smaller 
the inferred Mh/L becomes. 

In order to lift this degeneracy, therefore, it is essen- 
tial to detect the deviation of N{> L) from a power law. 
Since the cumulative mass function in 5.7 2 8.8 begins 
to decline exponentially with mass at a few times lO'^^ Mq 
at 2; ~ 6 to 10^^ Mq at 2 ~ 9 (see Fig. [TJ, an accurate 
determination of the bright end of the luminosity function 
at LLyc > W*^ erg s~^ x (100 LLya/Mh) at 2 ~ 9 to 
10*'^ erg s~^ X (100 LLya/Mh) at 2 ~ 6 would be required 
to lift the degeneracy. (Note that Mh/L is always quoted 
in units of MqLq, where Lq = 3.8 x 10^^ erg s"^) The 
bright end of the luminosity function is not constrained very 
well by the existing surveys (see Fig. [5)). One would need a 
larger survey area for a better determination of the bright 
end of the luminosity function. (See section 14.11 for further 
discussion on constraining the bright end of the luminosity 
function.) 

We have chosen to work with e = 1, which is allowe d 
by the existing data l|Diikstra et al.ll2006l : lstark et al.ll2007l '). 
As we show in this paper, this assumption does provide rea- 
sonable and useful results. We then use the degeneracy line, 
{Mh/L)e~^^^ = constant, to incorporate the effect of e into 
the inferred constraints on properties of Lya emitters. 

A further improvement to the model can be made by 
taking into account the fact that a relation between the lu- 
minosity and halo mass is not unique, but has some disper- 
sion. One may include this by using a conditional proba- 
bility of luminosity given the halo mass, P{L\Mh), as 7V(> 
L) = V{z) f dMh P{L\Mh) dn/d Mh with P(L\Mh) given 
by, e .g., a log-normal distribution (|Coorav fc Milosavlievid 
120051 ). Once again, the current data cannot constrain the ex- 
tra parameters characterizing P{L\Mh), except for its first- 
order moment, a mass-to-luminosity ratio. In our simpli- 
fied approach we take it to be a delta function, P{L\Mh) = 
5'^[Mh - {Mh/L)L], which gives equation 

In thi s paper we shall use the Sheth-Tormen form ula for 
dn/dMh (Sh eth et al.]|200ll : ISheth fc Tormen"2002') . Since 
the Press-Schechter mass function (Press & Schcchtcr 197|) 
tends to underestimate the number of haloes in the high 
mass range, the mass-to-light ratio inferred from the Press- 
Schechter mass function would be smaller than that from 
the Sheth-Tormen formula. We have found that the mass- 
to-light inferred from the Press-Schechter mass function is 
smaller by a factor of two. 

The volume, V{z), is found by multiplying the comov- 
ing volume element, d^V/dzdQ, by the depth of the survey 
(found by integrating over redshift) and the survey area on 
the sky. For narrow-band surveys, the redshift integral can 
be approximated as Az. The comoving volume element is 
given by 



dW _ cdl{z) I 1 



dzd^ l + z \ Hiz){l + z) 



where di is the proper luminosity distance. 



(2) 
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Figure 1. Mass function. The dotted, dashed, dot-dashed and 
solid hnes show z = 5.7, 6.55, 7.025, and 8.76, respectively. These 
are the redshifts of narrow-band surveys that we consider in this 
paper (see Table[TJ. ( Upper Panel) The Sheth-Tormen mass func- 
tion as a function of halo mass, dn/dM/^, in units of comov- 
ing Mpc"'' '^0^ {Lower Panel) The cumulative mass function, 

r°° dMfi dn/dMh, in units of comoving Mpc~^. The cumulative 
luminosity function can be calculated by stretching the horizon- 
tal axis of this figure by dividing Mf^ by a suitable factor that 
converts the mass to luminosity: a mass-to-light ratio, Mf^/L. 



2.2 Basic formalism 

From equation ([1} we can derive the number of 
observed above a certain flux density as: 



dz 



dFdh'^ 



d^V 
dzdQ 



dF 



dn dMh 



A2 



d^V 
dzd^l 



dMh dF 

it 

dn dMh 



Fur, 



dMh dF 



■&{Mh - M^,r.{z))dF. (3) 



where Mh is the lialo mass. 

Not all dark matter haloes will be forming stars - 
only haloes with a mass above some critical minimum 
mass {Mmin)- This is represented by the function ■d{Mh — 
Mmin{z)), which is zero if the halo mass is smaller than 
Mmin and unity if it is larger than or equal to Mmin- The 
minimum mass is only theoretically known, and we use the 
virial mass of a 10, 000 K halo, Mmin = 0.94 x 10* Mq[{1 + 



z)/W]-'^^^. However, given the current sensitivity of tele- 
scopes, it is unlikely that a halo of mass Mmin will be bright 
enough to be seen, unless the mass to light ratio of galax- 
ies is unusually small. Therefore, Mmin is irrelevant to our 
analysis presented in this paper and our conclusion is inde- 
pendent of the actual value of Mmin- 

In deriving equation we have made an assump- 
tion that each dark matter halo above Mmin hosts one 
galaxy. This is a valid assumption at high redshifts, as mas- 
sive haloes such as groups (M > 10^^ Mq) and clusters 
(M > 10^* Mq) of galaxies hosting multiple galaxies are 
extremely rare (see Fig. [1} . If we were to assume a field size 
equal to the largest survey area discussed in this paper (that 
of the LALA survey, 1296 arcmin^), and the widest redshift 
range (that of the Subaru survey at z~ 7), the number of 
haloes above a mass of 10^^ Mq can be found by using equa- 
tion [T] At a redshift of 1, there would be 32 haloes larger 
than this mass in the field. At higher redshifts (z = 3,5, and 
7), there would be less than one such massive halo in the 
field (0.98, 1x10"^, and 7.7 x 10"*, respectively.) Thus at 
high redshifts, it is safe to assume that there arc no groups or 
clusters observed. (Larger surveys (fOuchi ct al. 2005) have 
observed protoclusters) . This property makes it possible to 
model the luminosity function of high-z galaxies without 
complications ari sing from galaxy for mation processes. Some 
observations (i.e. lOuchi et al.l (|2004h ) show that some dark 
matter haloes host more than one Lya emitter. However, 
since galaxy occupation number at high redshifts is not well 
known, we will just assume one galaxy per halo. 

The most important uncertainty in our model is that 
not all galaxies are seen as Lya emitters. Some galaxies do 
not produce as many Lya photons as the others do because 
of dust extinction in galaxies themselves and scattering in 
the IGM. Therefore we could have assumed that there is less 
than one Lya emitter per halo; however, this effect can be 
modeled effectively by introducing a Lya survival fraction, 
Uesc- This parameter quantifies the fraction of Lya photons 
that escaped from a halo and the IGM. Therefore, as aesc 
increases, a galaxy is seen as more luminous intrinsically. 

In summary, we model the galaxy number counts by 
placing one galaxy per halo with only a fraction of photons 
escaping from galaxies and the IGM. 



2.3 Mass-to- "observed light" ratio 

Our formulation is now reliant on how we relate the flux 
density of a galaxy to its mass. 

The flux density of galaxies depends on two things; the 
luminosity distance to galaxies, and a mass-to- "observed 
light" ratio, which will relate the total mass of the halo 
(including dark matter) to the luminosity that is actually 
observed. The flux density of a galaxy observed by a certain 
instrument is found by: 



F : 



Lband/Mh 
ATld\{z)l\Vob, 



■Mh, 



(4) 



where Ai'obs is the bandwidth of the instrument (which we 
have assumed to have an ideal rectangular bandpass) and 



Lband ~ / dvLiy, 



(5) 



4 Elizabeth R. Fernandez and Eiichiro Komatsu 



is the observed luminosity within a certain bandwidth of the 
instrument, the rest-frame luminosity per unit rest-frame 
frequency, and 1/2 and ui the frequency limits of the survey. 

We assume that the mass to light ratio is independent of 
mass, dMh/dF = Mh/F. This approximation is well justi- 
fied, as the current surveys are probing a limited mass range. 
Total and stellar masses are still unknown for Lya emit- 
ters. Studie s of the stellar masses using SEP fitting have 
just begun (jMobasher et al ] I2OO5I : iNilsson et ai]|200'iD . In 
the future when the observations of Lya emitters can cover 
a wide mass range, one may use a parametrized model, e.g., 
, to improve fits. For the present purpose additional 
parameters are unnecessary. 

Using equation Q, equation Q can be rewritten with 
respect to Mn/Lband- 



Flir, 



dFdO. 



dF 



, ,2 , dV . . Mh 

= ATTdL{z)-z——/\zAVobs-f 

dzdil Ltanc 

Xi){Mh - Mminiz)). 



Flin 



dn{Mh{F)) 
dMn 



dF 
(6) 

To evaluate dn/dM for a given F, we use equation @ to 
convert F to Mh- Once again, Mh/Lband is independent of 
Mh, and Mh is almost always greater than Mmin, and thus 
almost always 'd{Mh — Mminiz)) ~ 1. The only unknown 
quantity in this equation is Mh/Lband, while Auobs, Fiimit, 
and Az are given by the survey properties. In other words, 
Mh/Lband is the parameter that should be measured from 
the observational data directly. 



3 NARROW-BAND LYa SURVEYS: 
OBSERVATIONS 

Several telescopes are now powerful enough to attempt to lo- 
cate galaxies at 2 > 6. One effective method of locating high 
redshift galaxies is to use narrow-band filters to detect Lya 
emission from a small range of redshifts. Lya emitters tend 
to have large equivalent widths, with their line intensities 
significantly higher than the continuum emission. In order 
to test that the galaxy is indeed a high redshift galaxy, it 
must be undetected at optical wavelengths, and follow-up 
spectroscopy may be employed if possible. At times, when 
no continuum is visible, the asymmetric profile of the Lya 
line can be used to identify a Lya emitter. 

A large number of such narrow-band Lya searches 
have been carried out on blank fields at z = 5.7 



jRhoads & Malhotral I2OOI ; Rhoads et all 


20031: Maier et all 


2003 


: Hu et al. 20041; Ouchi et al 


! 


20051: 


Aiiki et al 


20031. 


2004 


. I2OO6I: IShimasaku et al.ll200e 




Muravama et al 


200'(tl. 


z 


= 6.5 jKodaira et al.l |200^ 




jTaniguchi et al. 


I2OO5I: 



200 



and 



8.8 |Willis fc Courbinll2005l : IWiUis et all booilcubv et al.l 
l2007l ). lKurket all (|20oi r used an alternative technique, a 
slitless-grism spectroscopy survey, which has also yielded 
a successful detection of a Lya emitter at z = 6.5. 
iMartin fc Saw icki (2004) performed a multi-slit windows 
search at z = 5.7. 

In this paper we use the observational data from six nar- 
row band Lya searches: two from the Infrared Spectrometer 
and Array Camera (ISAAC) on the Very Large Telescope 



(VLT), three from Subaru, and one from the Mayall Tele- 
scope at Kitt Peak. The basic parameters of these surveys 
are summarized in Table [T] 

Two surveys at z = 5.7 and 2 — 6.56 were taken 
at the Subaru telescope using the Suprime-Cam. The sur- 
veys covered 725 and 875 arcmin^, respectively. These sur- 
veys had a limiting fiux density of 145 nj. Candidates were 
found down to a limiting magnitude of 26.0 (at a detection 
limit of 3(T for z=5.7 and 5a for z=6.56) using a 2" aper- 
ture. At 2 = 5.7, there were 89 Lya emitter candidates, 
63 were followed up by s pectroscopy, and 34 Lya emitters 
have been confirmed (iiShimasaku et aill2006f ). At 2 = 6.56, 
there were 58 Lya emitter candidates, 53 were followed 
up by spectroscopy, and 17 Lya emitters have bee n con- 
firmed l|Taniguchi et al.ll2005l : iKashikawa et al.ll2006l ). From 
these detections, they were able to fit a Schechter function, 
(j)(L)dL — (p* (L/L*)" exp{—L/ Li,)dL/L* , with parameters 
given in Table 12] 

The Large Area Lyman Alpha (LA LA) survey searche d 
for galaxies at a redshift of around 6.55 l|Rhoads et al.|[2004l ). 
This survey was conducted on the 4-m Mayall Telescope at 
Kitt Peak. A narrow-band 1296 arcmin^ image was taken of 
the Bootes field. The limiting flux density was 700 nJ. Three 
candidates were located down to a limiting magnitude of 
24.3 {5a) with a 1".02 aperture - one of which was confirmed 
spectroscopically as a Lya emitter. 

Another survey using the Subaru telescope covered an 
area of 876 arcmin'^ at 2 = 7.025. They used the narrow 
band filter NB973 on the Subaru Suprime-Cam and had a 
limiting flux density of 398 nJ. There were 2 Lya emitter 
candidates down to a limiting magnitude of 24.9 {5a) with a 
2" aperture - one of which w as confirmed sp ectroscopically 
as a Lya emitter at 2 = 6.96 (|lve et al ] l2006t ). 

ZEN, which stands for 2 equals nine, is a narrow J-band 
mission using the ISAAC on the VLT. Its central redshift is 
8.76 and it covers a field size of 4 arcmin^ down to a limiting 
flux density of 302 nJ. They searched for galaxies in the 
Hubble Deep Field South that displayed an excess in the 
narrow band in comparison to the J-band {Js — NB ^ 0.3) 
and that were undetected in the optical. No galaxies were 
found do wn to a limiting magnitude of 25.2 {5a) wit h a 0".7 
ape rture JWiUis fc Courb"inll2005l : IWillis et al.ll2006l ). 

ICubv et al.l l|2007^ did a followup narrow-band search, 
using the ISAAC at the VLT, with a larger held of view 
(hereafter referred to as the ISAAC ext). They imaged seven 
flelds in the Chandra Deep Field South that totaled 31 
arcmin^ down to a flux density of 1740 nJ. They also de- 
tected no galaxies within the flelds down to a limiting mag- 
nitude of 23.3 {5a) with a 1" aperture. 



4 PROPERTIES OF LYa EMITTERS 

4.1 Extracting a mass-to- "observed light" ratio 

Can we infer anything about properties of the Lya emitters 
discovered in these narrow-band surveys? Is the lack of de- 
tections in the flelds above 2 > 7 expected, or should we 
expect many more galaxies? 

In equation (|6|, the only free parameter was the mass- 
to- "observed light" ratio, Mh/Lband. Therefore, we vary 
Mh/Lband to give a model that is consistent with observa- 
tions. The number of galaxies drops as Mh/Lband increases. 
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Figure 2. The observed luminosity function of Lya emitters constrains their mass-to- "observed light" ratio. Each panel shows the 
cumulative number of sources detected in each field above a certain flux density, N{> F). (The flux density limits of each survey are 
indicated by the vertical lines with right arrows.) The bottom and top axis show the measured flux density and luminosity (in erg s~^), 
respectively, while the left and right show the number of sources per field and per comoving Mpc'^ , respectively. The mass-to- "observed 
light" ratio, M^/Li,^„^, is fit to each of the luminosity functions derived from various narrow-band surveys. Here, Lfcand refers to the light 
that falls within the band of instruments, which is mostly the Lya line. Curves in each panel show predictions with various Mf^/ Lf,and- 
{Upper Left) The Subaru field at z = 5.7 jShimasaku et al.|[20o3) . The stars, triangles and diamonds show their Schechter-fit to the 
luminosity function with a = —2, —1.5 and —1, respectively. The scatter between symbols shows uncertainty, and they diverge mostly 
below the flux density li mit, as expected. T he solid circles are the data, with the error bars showing Poisson error. ( Upper Right) The 
LALA field a.t z = 6.55 l lRhoads et al.ll2004l. One Lya emitter was found, and th e horizontal line shows A'^(> _F) = 1 per field. {Middle 
Left) The Subaru field at z = 6.56 l|Taniguchi et al.ll2005l: iKashikawa et al.ll20od) . The meaning of symbols is the same as in the upper 
left panel. Data is from both the photometric and spectroscop ic samples. The circles show data corrected for detection completeness. 
{^Middle Right) The Subaru field a.t z = 7. 025 dive et al.ll2006ll . One Lya emitter was found. {Bottom Left) The ZEN field at z = 8.76 
Iwillis fc Courbinll2005l : IWillis et al.ll2006l) . No sources were found , and the horizontal line with lower arrows shows N{> F) < 1 per 
field. {Bottom Right) The ISAAC ext field at 2 = 8.76 jCubv et al.ii2007i) . No sources were found. 
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Tanieuchi et alj \200^ Y. 
Cubv et alj l|2007h ''. 



Table 1. Survey par a meters taken fro m IShimasaku ct al.' ("ioOff)", 'Rhoads et al.' ('2004')', 
iKashikawa et al.l ll2006h''. live et al.1 1|2006|)'=. Willis & Courbin (2005)-'', Willis et al. (2006)!?, and 



Name of survey 


Telescope 


Central A (A) 


Bandwidth (A) 


Central z 


Redshift Range 


Area (arcmin^) 


Ref 


Subaru Deep Field 


Subaru 


8150 


120 


5.7 


5.64 - 5.76 


725 


a 


LALA 


Mayall 


9182 


84 


6.55 


6.516 - 6.586 


1296 


b 


Subaru Deep Field 


Subaru 


9196 


132 


6.56 


6.508 - 6.617 


876 


c, d 


Subaru Deep Field 


Subaru 


9755 


200 


7.025 


6.94 - 7.11 


876 


e 


ZEN 


VLT 


11900 


89.5 


8.76 


8.725 - 8.798 


4 


f, g 


ISAAC ext 


VLT 


11900 


89.5 


8.76 


8.725 - 8.798 


31 


h 



Table 2. Best-fitting Schechter parameters from IKashikawa et al] (l2006h . 
These data are plotted in Figure [5] 



Redshift 


a 


^oZwiL*/h^Q erg s ^) 


\og^o{r /h^O Mpc-3) 


6.5 


-2.0 


42.74 


-3.14 




-1.5 


42.60 


-2.88 




-1.0 


42.48 


-2.74 


5.7 


-2.0 


43.30 


-3.96 




-1.5 


43.04 


-3.44 




-1.0 


42.84 


-3.14 



(See Fig. O) As the mass-to-light ratio increases, the star 
formation is spread out over a longer period of time. There- 
fore, the galaxies are dimmer because less stars are shining 
at any given time. 

We find that the Subaru data at z = 5.7 and z = 6.56 
are consistent with no evolution of properties of Lya emit- 
ters or the IGM opacity. The evolution in the number density 
of Lya emitters can be explained solely by the evolution of 
the halo mass function. These poin t s have been made a l- 
rcady bv'Malhotra fc Rhoadsl (|2()04h : jPiikstra et all (|20o3); 
[ivIcQuinn ct al. ( 2006). At 2 = 5.7 we find that Mn/Lband = 
95 — 120 fits the Subaru data from lShimasaku et al. I (|2006l ). 
with lower values favored near the flux density limit, where 
it is harder to correct for sample completeness. At z = 6.56 
a slightly (20%) smaller val ue, Mh/Li,and = 85 — 100, fits the 
Subaru data from jKashikawa et al.l (|2006i ). (These are fit for 
the values of Ltand below 10''^ erg s~^ to the flux density 
limit, for most of the observational data fall between these 
limits.) 

The high value of Mh /Ltand = 95 - 120 at z = 5.7 
seems to fall outside the Poisson error at z = 6.56 near the 
flux density limit, although at higher fluxes its statistical 
significance is more questionable. If we take this ~ 20% de- 
crease in Mh /Lband seriously, an interesting conclusion may 
be drawn. First of all, the decrease in Mh./Lband from z = 5.7 
to z — 6.56 is qualitatively inconsistent with the evolution 
of neutral fraction in the IGM. If the IGM was more neutral 
(i.e., less ionized) in the past, we should observe the increase 
in Mh/Lband at higher z. The evolution in M^/Lband may 
be even more significant than it looks now, once the redshift 
effect is taken into account. Since the survey at z — 6.56 
collects less photons than that at z = 5.7 for a given band- 
width of the instrument, one must take into account the 
bandwidth properly before making a quantitative compari- 
son between M^/ Lband from two different redshifts. We shall 
perform this analysis more carefully in section 14.41 



Our finding may suggest that (i) Lya emitters at z = 
6.56 are brighter intrinsically than those at z = 5.7, or (ii) 
the intrinsic luminosity is the same, but more Lya photons 
escaped from galaxies at z — 6.56 than from z = 5.7. (The 
absorption in the IGM was kept the same.) The possibility 
(ii) is quite plausible, if dust content of galaxies at z = 6.56 
is less than that at z — 5.7. How much less requires a more 
careful analysis, which we shall give in section IT4l 

In order to put better constraints on the bright end 
of the luminosity function, a larger survey is needed. The 
brightest Lya emitter detected at z=5.7 had a narrow band 
magnitude of 23.41 and the brightest at z=6.56 had a narrow 
band magnitude of 24.13. To detect brighter Lya emitters, 
a larger survey would be needed that could find the rare, 
high density peaks of the mass function. In order to detect 
10 galaxies above a magnitude of 23.41 at z=5.7, the sur- 
vey would need to be at least 6.34 degree'^ if we assume 
Mh/ Lband of 120, and at least 3.13 degree^ if we assume a 
Mh/ Lband of 95. For the Subaru field at z=6.56, to detect 10 
galaxies with a narrow band magnitude of 24.13 or higher, 
the area of the survey would have to be 3.6 degree^ if we 
assume a Mh/ Lband of 100 and 2.2 degree^ if we assume a 

Mh/ Lband of 85. 

The LALA data l|Rhoads et al.ll2003 ) also probe a very 
similar redshift, z = 6.55. Since only one galaxy was found 
from LALA, the Poisson error is large. Nevertheless, the 
LALA data give us an important cross-check of the results 
obtained from the Subaru field at the same redshift. We 
find that Mh/ Lband ~ 160 explains LALA's detection of one 
galaxy at z = 6.55. When we compare the LALA and Subaru 
counts, we must take into account the different bandwidths 
of these surveys. The LALA's bandwidth is about 60% nar- 
rower than Subaru's (Table [T]), and thus the constraint from 
the LALA data would correspond to Mh/ Lband ~ 100 for 
the Subaru data. We thus conclude that the constraints from 
the Subaru and LALA fields at z — 6.55 are comfortably 



Mass-to-light Ratio of Lya Emitters 7 



consistent with each other. A more thorough comparison 
will be given i n section 14.31 

At z = 7.0 live et al.l ( 2006h discovered one Lya emitter 
in the Subaru field that was conformed spectroscopically. 
While the Poisson error is large, we find that Mh/Ltand ~ 
100 explains Subaru's detection of one galaxy ai z — 7.0. 
This number is remarkably similar to what we have found 
from the Subaru fields at z = 5.7 and 6.56 as well as from 
the LALA field at z ^ 6.55. 

At z = 8.76 the searches performed in the ISAAC/VLT 
fields yielded null results. We therefore place lower limits 
to Mh/Lband > 7 and 3 from the ZEN and ISAAC ext, re- 
spectively. The weaker constraint from the latter is due to a 
brighter flux density limit. These lower limits are consistent 
with properties of Lya emitters as constrained by the other 
searches at z ^ 7. 

In order to improve upon these results, how large of a 
survey would be needed? In order for 10 galaxies above the 
flux density limit to be seen in the Subaru (z=7.025) and the 
LALA flelds, the survey area would need to be increased to 
2.5 degree^ and 3.4 degree'^ respectively. If we assume that 
the mass to light ratio and aeac is the same at z=8.76 than 
it is at z=7.025, the survey area of the ZEN and ISAAC 
ext fields would need to be increased to 24.8 and 1.69 x 10^ 
degree^ respectively. The area needed for the ISAAC ext 
field is larger than the entire sky, so it becomes apparent 
that an adjustment to the bandwidth or fiux density detec- 
tion limit is necessary to make finding a Lya emitter more 
feasible. 

In summary, properties of Lya emitters and the IGM 
opacity have not evolved very much between z = 5.7 and 
7. The lack of detection at z — 8.76 is also consistent with 
no evolution, although it does not provide a significant con- 
straint yet. 



4.2 Finding a mass-to- "bolometric light" ratio 

Our analysis so far has been relatively model-independent. 
We have extracted the only free parameter, Mh/Ltand, from 
various narrow-band searches of Lya emitters. 

Here, Mh/Lband only describes the light observed over 
the narrow band (the luminosity within the bandwidth, 
Lband, not the bolometric luminosity, Lboi)- To proceed fur- 
ther and understand physical properties of Lya emitters 
better, however, we must relate Mh/Lband to the mass-to- 
"bolometric light" ratio, Mh/Lboi, taking into account stel- 
lar populations as well as differences in the bandwidths. 

In order to get the actual mass to light ratio, the spec- 
tra of a stellar population of galaxies must be modeled and 
integrated first over all frequencies and then compared to 
the light that is observed in the narrow band. The fraction 
of Lyof photons that survived, Qesc, also needs to be taken 
into account. 

As a result, each data-set yields a constraint on 
(Mh/Lboi) function of assumed stellar populations. 

When the duty cycle is less than unity, the constraint 
should be interpreted as {Mh/ Lboi){ciEsct^^'')~^ , where 7 ~ 
2 is a local slope of the cumulative luminosity function, A'^(> 
L) oc L to which the current data are sensitive. 

The spectra of a population of stars depend on the stel- 
lar mass spectrum and metallicity of stars. We shall use a va- 



riety of mass fu nctions paired with metallicities: (a) Salpeter 
(|Salpeteij|l9"55l ): 



rr \ -2.35 

j(rn) oc m , 

(b) Larson (|Larsonlll99l l: 



/(m) (X m 



"(1 + -)" 



(7) 



(8) 



which matches Salpeter's in the limit of nic 0. One can 
explore a variety of models by changing one parameter, rric. 
(c) A top-heavy spectrum: 

100 < m < 500 Mq 



/(m) (X 



0, 



otherwise, 



(9) 



which might be possible for the primordial metal-free stars 
(|Bromm fc Larsonll2004l '). (Note that m/(m) is fiat for 100 < 
m < 500 Mq.) The normalizations are given by 



dm f(m) 



(10) 



with mi and rrii being the mass limits that the mass function 
is integrated over. We also consider a delta function mass 
spectrum, with populations consisting of only 300 Mq stars. 

The synthetic spectrum emerging from a galaxy with a 
given population of stars is the result of a variety of radiation 
processes. Some of the light from the star is converted by 
the nebula into the Lya line, free-free, free-bound, and two- 
photon emission. 

We use analytical formulae for these spectra given in 
section 2 of iFernandez fc Komatsd ||2006|). paired with a 
line profile o f Lya emission from Loeb fc Rvbickil |l99^; 
ISantos et ahl (|2002). These formulae are fully analytic and 
therefore it is easy to adjust various parameters. We have 
checked that our predicted luminosities agree well with those 
produced from the numerical code, cloudy. 

Care must be taken when one computes the profile of 
Lya line. The width of the line profile we adopt here is likely 
too broad, as this profile assumes that the IGM around a 
source is completely neutral. The other extreme delta 
function profile at 10.2 eV, increases the inferred Mh/Lboi 
by a factor of at most a few. We study this issue further in 
Appendix 1X1 

We explore three metallicities: what we refer to as Pop- 
ulation III {Z = 0), Population II {Z = 1/50 Zq), and Pop- 
ulation I {Z — 1 Zq). For convenience we have fit the theo- 
retical stellar data for bolometric luminosity, stellar temper- 
ature, lifetime, and number of ionizing photons. The fitting 
functions, given in Table [3] and plotted in Fig. |3l have been 
obtained by fitting stellar models provided by the papers 
given in the third column of Table O 

In Fig.|4l we show the luminosity integrated over a mass 
spectrum divided by the average stellar mass: 



(11) 



for each component i, which includes stellar black-body, 
Lya, free-free, free-bound, and two-photon. 

Now we are in a position to calculate a conversion factor 
from Lband to Lboi- 

Lband 



Lboi ir "^""T^.n^ L^.^fi^)dm 



(12) 
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Table 3. Fitting functions for the number of hydrogen ionizing photons per second, Q{H), stellar temper- 
ature of the star, T^ff, bolometric luminosity of th e star, L*^^, and stel lar lifetime, t<,, for varying m etallic- 
ities. These wer e obtained from stellar models from Marigo ct alj l|200ll ')°.[Lei eune fc Schaereij l|200ll ')''. and 
ISchaered l l2002h '^. or the fitting functions from iSchaerer, 120021') ''. Note that y = log(Mt/M0), and log is a 
logarithm of base 10. 



Z = 


Fitting Function 


Ref. 


log[Q{H)/s-i] 


43.61 + 4.90y - O.SSj/^ for 9 ^ M. ^ 5OOM0 


d 




39.29 + 8.551/ for 5 ^ A/, < 9Mq 


d 




for Af, < 5 




log[T,^^/K] 


3.64 + 1.50?/ - 0.556?/2 + omOy^ for M, ^ IOMq 


c 




3.87 + 0.937?/ - 0.156?/2 for M, < IOMq 


a 




0.457 + 3.901/ - 0.5301/2 j^^, ^ IOMq 


c 




0.219 + 4.511/ - 0.9231/2 for A/, < IOMq 


a 


log [t* /year] 


9.79 - 3.76?/ + 1.411/2 _ Q igQyS for M* ^ IOMq 


c 




9.80 - 3.52?/ + 1.181/2 - 0.164?/=' for Af* < IOMq 


a 


Z = 1/50 Zq 


log[Q(H)/s-i] 


27.80 + 30.68y - 14.80?/2 + 2.50i/3 for Af* > 5 


d 




for Af, < 5 




log[T,^^/K] 


3.86 + 0.898?/ - 0.299?/2 + 0.0399?/3 


b 


^osILIoi/Lq] 


0.175 + 4.261/ - 0.6521/2 


b 


log [r* /year] 


9.82- 3.13?/ + 0.743y2 


b 


Z = 1Zq 


log[Q(J?)/s-l] 


27.89 + 27.751/ - 11.87?/2 + 1.73i/3 for Af, ^ 5 


d 




for Af, < 5 




log[T,^^/K] 


3.74 + 0.826?/ -0.1661/2 


b 


logl^L;/-^©] 


-0.148 + 4.70?/ -0.781?/2 


b 


log [t* /year] 


10.08 - 3.47?/ + 0.7741/2 o.0327i/3 


b 
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Figure 3. Fitting functions for the theoretical stellar data for the stellar bolometric luminosity, i^^p lifetime, 
r, , the number of hydrogen ionizing photons per second, Q{H), and the stellar temperature, T^ff, for varying 
metallicities: Z = (Population III), Z = 1/50 Zq (Population II), and Z = 1 Zq (Population I). The 
fitting formulae are given in Table [3] In addition, the nuclear burning efficiency (the total radiation energy 
produced over the stellar lifetime divided by the rest mass energy, Af, c2), and the total number of ionizing 
photons over the star's lifetime are plotted. 



Mass-to-light Ratio of Lya Emitters 



Free-free 

Free-bound 

Stellar 

Lyman olpha 

Two-photon 

Totol 



Luminosity for the Heovy Moss Function, Pop 




Lunnlnosity for the Lorson Moss Function, Pop III, = 50 



hy [ev] 

Luminosity for the Lorson Moss Function, Pop III, = 10 




hi/ [eV] 

Luminosity for the Lorson Moss Function, Pop II, m^ = 50 



. EW = 577 



hi/ [eV] 

Luminosity for the Solpeter Moss Function, Pop I 




hv [eV] 

Luminosity for the Lorson Moss Function, Pop I, = 10 




hy [ev] 



EW = 1652 



j: 



hi/ [eV] 

Luminosity for the Solpeter Moss Function, Pop I 




hi/ [eV] 

Luminosity for the Lorson Moss Function, Pop II, m^ = 10 




hi/ [eV] 

Luminosity for the Lorson Moss Function, Pop I, m,. = 50 



EW = 262 



J 



? 21 ^ ' -..-<f^ 



hi/ [eV] 

Luminosity for the Solpeter Moss Function, Popi 




hi/ [ev] 



Figure 4. Rest-frame spectra of galaxies with various populations of stars integrated over a mass spectrum. 



The vertical axis is in units of erg s ^ Hz 



M 







The solid lines show the total spectra, while the dotted, 



short dashed, dot-dashed, dot-dot-dot-dashed, and long dashed lines show the free -free, free-bound, stellar , 
Lyg, and two-phot on emission, respectively. We have adopted a Lya profile from iLoeb fc Rvbickil l|l999l '): 
ISantos et al.l 1I2OO2I) . (For further discussion of the profile, see AppendixjX]) We show the luminosity averaged 
over various mass spectra, given in section |4l divided by the mean stellar mass. The spectra are computed 
for a galaxy at 2: = 7.025, although the redshift affects the shape of the Lya line profile only. The equivalent 
width (EW) of the Lya line before extinction or scattering is also given. The EW has been computed from 
EW = (Total flux in Lya)/(Continuum flux at 1216 A). 
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Table 4. Ratio of luminosity observed within the bandwidths of Subaru (L5ub,z=5.7i Shimasaku et al. 2006; 1/3116,2=6.561 
Taniguchi et al. 2005; Kashikawa et al. 2006; Lsub,z^7, lye et al. 2006), LALA (Llala, Rhoads et al. 2004), and ZEN 
and ISAAC ext (Ljsaac, Willis & Courbin 2005; Willis et al. 2006; Cuby et al. 2007) in comparison to the bolometric 
luminosity, for a variety of mass spectra and metallicities, in the absence of absorption or extinction of Lya. The ratio 

is nearly independent of stellar mass spectra, while it drops as the metallicity of stars increases. 



Metallicity {Zq) 


f(m) 


mi,m2(MQ) 










LiSAAC 








Lbol 







3OOM0(5-function 




0.312 


0.209 


0.263 


0.297 


0.120 





Heavy 


100, 500 


0.311 


0.209 


0.262 


0.297 


0.120 





Larson, rric = 50Mq 


0.8, 150 


0.313 


0.210 


0.264 


0.299 


0.120 





Larson, rric = IOMq 


0.8, 150 


0.305 


0.205 


0.257 


0.291 


0.117 





Salpeter 


0.8, 150 


0.293 


0.197 


0.247 


0.280 


0.113 


1/50 


Larson, rric = 50Mq 


0.8, 150 


0.219 


0.146 


0.185 


0.210 


0.0843 


1/50 


Larson, rric = IOMq 


0.8, 150 


0.209 


0.139 


0.176 


0.201 


0.0806 


1/50 


Salpeter 


0.8, 150 


0.198 


0.131 


0.167 


0.190 


0.0761 


1 


Larson, nic = 5OM0 


0.8, 120 


0.143 


0.0945 


0.121 


0.138 


0.0555 


1 


Larson, nic = IOMq 


0.8, 120 


0.134 


0.0881 


0.113 


0.130 


0.0519 


1 


Salpeter 


0.8, 120 


0.124 


0.0813 


0.105 


0.120 


0.0480 



where Vmax = !^2,o6s(l + z) and Vmin = vi.obsi'^ + z) axe 
the limiting frequencies of instruments in the rest frame of 
the galaxies (again assuming a rectangular bandpass), Lv,i 
is the luminosity of each component, and m2 and mi are the 
mass limits of the stellar mass spectrum, /(m). 

In Table |4] we show Lband/ Lboi for the various surveys. 
We find that Ltand/ Ltoi is fairly constant over different mass 
spectra, but depends mainly on metallicity. Metal-free stars 
sustain higher temperatures as they undergo nuclear burn- 
ing through the p-p chain. Because of this, their stellar spec- 
trum is harder than stars with metals (see panels a-d in 
Fig. |4)| . Therefore they emit more ionizing photons that can 
be converted by the surrounding nebula into the Lya line. As 
the metallicity increases (panels e-j), the stellar temperature 
decreases, the stellar spectrum softens, the ionizing photon 
flux decreases, and thus the Lya line is depleted for a given 
stellar mass. As a result, one obtains a lower Ltand/ Ltoi for 
a higher metallicity. 

Now, calculating the mass-to-bolometric light ratio of 
galaxies is simple: multiply the mass-to-observed light ratio 
{Mh / Lhand) by the ratio of observed to bolometric luminos- 
ity given in Table [J] 

This is, however, not the end of story. Not all Lya pho- 
tons would escape from galaxies due to dust extinction, or 
from the IGM due to resonant scattering. We are unable to 
distinguish these two effects; thus, we parametrize a com- 
bined effect by a single parameter, a^sc, a survival fraction 
of Lya photons. 

4.3 Results 

The outcome of our analysis is a mass-to-bolometric light 
ratio divided by a Lya survival fraction and the effect of 
duty cycle, {Mh/ Lboi){oLesc(^^'')~^ , where Lboi is the intrin- 
sic luminosity of galaxies before absorption or extinction of 
Lya photons. We tabulate this quantity inferred from var- 
ious narrow-band searches in Table [5] which is the main 
result of this paper. 

Since the effect of bandwidths has been taken into ac- 
count properly, these constraints can be compared with each 
other on equal footing. Let us analyze the results in Table 



[S] As noted in the previous section, the Lya line diminishes 
in strength and [Mu/ Lboi){ciesc(^^'')~^ drops as metallicity 
increases. A variation due to different stellar mass spectra 
is negligible. 

The values of {Mu/ Lboi){ciesc(^^~')~^ inferred from the 
current data processed through our simple model are rather 
reasonable: for all cases where at least one source is found 
per field, the inferred {Mh/ Lboi){oi<:ac(^^'')~^ falls between 9 
and 38, the low and high values being for the solar and zero 
metallicity, respectively. Consistency across redshifts (2 = 
5.7, 6.5, and 7.0) as well as across different observations is 
striking. 

We conclude from these results that the Lya emitters 
detected in these narrow-band surveys are either normal 
galaxy populations with Mh/Ltoi ~ 10 and having a fair 
fraction of Lya photons escape, aescf^''^ ^ 0.5 — 1, or star- 
burst galaxies with Mh/ Lboi ~ 0.1 — 1 and a smaller fraction 
of the Lya photons escaped from the galaxies themselves and 
the surrounding IGM, a^scf^^'^ ~ 0.01 - 0.05. Note that the 
degeneracy still allows for a possibility of having a signifi- 
cant survival fraction from these starburst populations, e.g., 
aesc ~ 0.5, \i e^/^ - 0.1 (or e ~ 0.01 and 7 ~ 2). 

It is clear that Oesc, e, and Mh/L^oi are completely de- 
generate: we can only constrain the product of these three 
properties, not the individual properties. 

The ratio of the Lya flux to the continuum fiux helps to 
lift this degeneracy partially. The observed equivalent width 
(EW) of Lya emitters is on the order of 100 A or larger (e.g., 
iKashikawa et al.ll2006l ). In Figure |3] we show the predicted 
EW, 

EW= Total Flux in Lya ^^^^ 

Continuum Flux at 1216 A 

Using the zero- age main sequence values for the luminosity 
(given in Table[3]), we obtain equivalent widths of 1300—1700 
for Z = Q, 480 - 580 for Z = 1/50 ^q, and 220 - 260 for 
Z = 1 Zq. Therefore, a low survival fraction, a^sc ~ 0.1, is 
required for low metallicity populations, while a high a^sc ~ 
0.5 is required for high metallicity ones, in order to fit the ob- 
served EW. However, as the age of stars within the galaxy in- 
crease, the fraction of ionizing photons to non-ionizing pho- 
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Table 5. The mass (total halo mass) to light (bolometric luminosity) ratio times l/(«esc e^/'''). The lumi- 
nosity refers to the intrinsic luminosity before absorption or extinction of Lya photons. For each metallicity 
of stellar populations (Z = 0, Z = 1/50 Zq, and Z = 1 Zq) a range of values represent a range of stellar 
mass spectra. For the observational data we used Mh/Ltand = 95 — 120 and 85 — 100 for the Subaru fields at 
2 = 5.7 and 6.56, respectively, whereas we used Mf^/ L^and ~ 160 and ~ 100 for the LALA field at 2 = 6.55 
and the Subaru field at 2 = 7.025. The latter values are much more uncertain than the former ones due to a 
large Poisson error, as only one Lya emitter was found in each of the latter fields. For the ZEN and ISAAC 
fields only lower limits are given, as no sources were found in these fields. Note that 7 ~ 2 for the surveys 
listed here. 



Field Redshift ^ ^-rj- (Z = 0) ^7- (Z = 1/50 Zq) 4^ ^7- (Z = 1 Zq) 



Subaru 


5.7 


28 - 38 


19 - 26 


12 - 17 


LALA 


6.55 


~ 32 - 34 


--21-23 


~ 13 - 15 


Subaru 


6.56 


21 - 26 


14 - 19 


8.9 - 12 


Subaru 


7.025 


~ 28 - 30 


~ 19 - 21 


~ 12 - 14 


ZEN 


8.76 


> 0.79 - 0.84 


> 0.53 - 0.59 


> 0.34 - 0.39 


ISAAC ext 


8.76 


> 0.34 - 0.36 


> 0.23 - 0.25 


> 0.14 - 0.17 



tons decrease, and thus less photons are converted into Lya 
photons. Therefore, the equivalent width may decrease with 
time - dependin g on the age of the galaxy and the rate of 



1 11-^^ II ^ 

star formation. (ICharlot fc Falllll993l: ILeitherer et aLlll999l: 

Kudritzki et~aLl I2OO0I : iMalhotra fc Rhoadd |2002| : ISchaerer I 
20031 ) 



Assuming the zero-age main sequence luminosity of the 
stars, there are two solutions left for e ~ 1: (i) Lya emit- 
ters at 2; 5.7 are normal populations with Z > 1/50 
and Qesc > 0.5, or (ii) they are starburst populations with 
Z < 1/50 and Oesc < 0.1. As stars age, the EW of the 
population will also decrease, allowing for larger values of 
Qesc for a. Z < 1/50 Zq population. For e < 1 other solu- 
tions are still allowed. 

Having the Lya line be diminished in flux by about 
an order of magnitude is not a surprising effect. Both the 
IGM and galaxies themselves are expected t o sca tter or 
absorb Lya photons efficiently. iDiikstra et al.l (|2007l ) claim 
that the asymmetry in Lya lines that has been seen in the 
current data already suggests that the IGM only transmit- 
ted 10 — 30 % of the Lyg flux. Several authors model th e ef- 
fect of dust jHansen k. Ohll2006l:IVerhamme et al.ll2006l) and 
neutral hydrogen ( Laursen fc Sommer-Larsen 2007h within 
galaxies on Lya photons, and find that even a small amount 
of dust can easily absorb Lya photons, and the resulting 
line profiles may be complex due to a structure in the distri- 
bution of dust and outflows. In addition, high opacity near 
the line centre decreases the flux at the centre of the line. 
However, if the medium is clumpy, the ratio of Lyg to con- 
tinuum photons might actually be in c reased (|Neufeldlll99ll : 
[Hansen fc Obll2006l ). iMcQuinn et al.l (|2006h also study the 
effect of neutral hydrogen on the Lya line using cosmolog- 
ical simulations. Once the universe is almost totally reion- 
ized, there will not be much suppression of the Lya line, but 
before then, Lya luminosity might be able to help probe the 
size of HII bubbles - the larger the bubble , the less suppres- 
sion of the Lya line l|Haiman fc Cenll2005l ). 

The physics of this problem is complex; however, our 
results are consistent with a depletion of the Lya ffux 
if indeed these Lya emitters are starburst galaxies with 
MulUoi - 0.1 - 1. 



4.4 Interesting features 

Is there any "anomaly"? Let us focus on the Subaru 
fields at z = 5.7 and 6.56, as these are the most ac- 
curate data-sets. We observe nearly 20-30% decrease in 
{Mh/Lboi){aesc^^^^)~^ from z = 5.7 to z = 6.56. We saw 
this trend in the preliminary analysis based upon Mh / Ltand 
in section |4T] After a more careful analysis we still observe 
the same trend. 

Although subtle, if this is indeed a real effect, what 
would be the implication? This effect cannot be explained 
by having a smaller Uesc (hence a larger opacity for Lya pho- 
tons) at higher 2. Therefore, it is inconsistent with neutral 
fraction in the IGM around sources being higher at higher 
z. On the contrary, one needs to have a larger a^sc - hence 
a smaller opacity for Ly a photons - at higher z, perhaps 
due to less dust content (|Haiman fc Spaanj[l999l ). An al- 
ternative possibility is that Mh/Lboi was lower in the past, 
i.e., the Lya emitters were intrinsically brighter at higher 
z, perhaps due to a more intense starburst. Such a burst 
would create a large HII bubble around the source, which 
also helps to increase a^sc by suppressing the IGM opacity. 
It therefore seems easy to explain the 20-30% decrease in 
{Mh/ Lboi){ceesc^^^^)~^ from z = 5.7 to 6 . 5. A s imilar trend 
has also been pointed out bv lStark et all |20o3). 

As we show in AppendixlX] the magnitude of this effect 
is reduced to 10-20% if we assume that a line profile of Lya 
photons is a delta function at 10.2 eV. 

Another interesting feature in Table [5] is that 
{Mh/ Ltoi){c(esc^^''^)~^ &t z — 6.56 for Z = agrees 
with that a.t z — 5.7 for Z = 1/50 Zq, and 

{Mh/Lboi)iaesce^^'^)~'^ at 2 = 6.56 for Z = 1/50 Zq agrees 
with that at 2 = 5.7 for Z = 1 Zq. While it seems a mere 
numerical coincidence, it might also be suggestive of the 
metallicity evolution in Lya emitters. 



5 COMPARISON WITH PREVIOUS WORK 

A halo mass function as a tool for calculating the lumi- 
nosity function of Lya emitters is not a new idea (e.g., 
iHaimari fc Spaans|[l999l : iHaiman et allbOOOl ). 
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Novelty of our approach is the use of the mass-to-hght 
ratio as a fundamental parameter, which has a few advan- 
tages. In this section we make this point clear by comparing 
our results with recent work on a similar subject. 

In addition, our analysis is new in that we have explored 
various assumptions about metallicity and stellar mass spec- 
tra of Lya emitters. 

5.1 iDiikstra et al.l (|2006l ;) 

iDiikstra et al.1 (120061 ) computed A''(> F) by integrating the 
halo mass function over mass above a certain flux density, 
F. In order to relate the host halo mass to the observed 
luminosity. La, they used 

^=0.128x^^4(15^, (14) 

which is their equation (2) in our notation. (Note that in 
our notation M/L is always measured in units of MqL^ .) 
Here, rj is the fraction of baryon mass converted into stars, 
tsys = ^thub is the duratiou of a starburst, thub is the Hubble 
time, and e is the duty cycle. 

In our approach Mh/La is the only free parameter, 
and the effect of e is included using the degeneracy line, 
{Mh/ La)e~^^'' — constant. Their approach was to divide 
Mh/La up further by introducing two free parameters, e 
and rjUcsc, and constrain these parameters simultaneously. 
However, it is difficult to extract more tha n Mh/Lg from 
the ob served luminosity function. Figure 1 of lDiikstra et al.l 
l|2006ll also shows that e and ryOesc are strongly degenerate. 
In our opinion the current data do not allow for two free 
parameters to be constrained well. In addition, the use of 
Mh/Lct as a parameter avoids the need to specify the dura- 
tion of a starburst or the fraction of baryon mass converted 
into stars. 

They found that Oesc at z = 5.7 inferred from 
Shimasaku et al.1 ||2006| ) and Qesc at z = 6.5 in ferred from 
Taniguchi et al.1 (|2005| ): iKashikawa et al.1 l|2006l ) are about 
the same, the ratio of the two being aesc.sr/aesc.es ~ 
0.8 — 1.5 for a prior on e of e = 0.5 — 0.03. We would find 
aesc,57/aesc,65 ~ 0.7 — 0.8 (scc Table [5]) , if we assumed that 
the intrinsic properties of Lya emitters did not change be- 
tween these redshifts. The other values of a^sc are permit- 
ted when we vary the intrinsic mass-to-light, Mh/Lboi, with 
{Mh/ Li,oi){aesct^^'^)~^ held fixed. This is essentially equiv- 
alent to their varying e along the degeneracy line. 

When an additional constraint from the luminosity 
function of UV continuum was included in the analy- 
sis, they found that the constraints shifted slightly to 
ciesc,57 / Q.<isc,&5 ~ 1-1 — 1.8. While we do not perform a joint 
analysis with the UV continuum luminosity function in this 
paper, we would expect a similar shift in the parameter con- 
straint. 

5.2 ISalvaterra &: Ferraral (|2006l ') 

ISalvaterra fc Ferraral l|2006l 'l used an stellar mass spectrum 

^ Their luminosity, La, is different from our Li,and, as they ig- 
nored the line profile, continuum, and bandwidth of instruments. 



of stars that is given by a delta function at m« = 300 Mq, 
and related Mh to Lband as 

Mh ^ 300 Mq l_ 

Lband XT""" di^E,'^"'." (300 Mq) i^??' 

which can be obtained from their equation (3), combined 
with our equation Q. Their Mh/ Lband is therefore equal 
to about 10 timefl the mass-to-light of a metal-free star of 
300 Mq. 

The lifetime of a starburst of their model galaxy is as 
short as the lifetime of stars, which is only 2 Myr. In other 
words, they assumed that these massive stars formed at once 
in a galaxy, so that the lifetime of starbursts was the shortest 
possible time, equal to the lifetime of the star. This creates 
a very short lived but extremely bright galaxy that could 
easily be detected with current observations. 

The mass-to-bolometric light ratio of their model galaxy 
was Mh/Lboi = 6.73 x 10"" and 1.35 x 10"^ for "H-cooling" 
and "H2-cooling" haloes, respectively, assuming all Lya pho- 
tons escaped. (They used r; = 0.8 and 0.4 for H-cooling and 
H2-cooling haloes, respectively.) These extreme values allow 
them to predict that there should be thousands of galaxies 
seen in the NICMOS Ultra Deep Field, where only three or 
fewer were actually detected, and 400 to 700 in the ZEN 
field, where no sources were detected. They also reported 
that almost all of the Spitzer counts should be attributed to 
galaxies above z ~ 8. 

Their conclusion is driven by their fixed value of 
Mh/Lijoi, which seems rather low. Our formulation, which 
treats Mh/Lf,and as a free parameter, allows for dimmer 
galaxies by spreading out the star formation over a much 
longer period than the stellar lifetime. This allows us to ob- 
tain results that are consistent with observations. Note that 
their using a delta-function mass spectrum is not the source 
of discrepancy. We can still fit the observations with a rea- 
sonable Mh/Lboi for the same mass spectrum. The source 
of discrepancy is their assumption about an instantaneous 
starburst in 2 Myr. 

They used these bright galaxies to fit the observed ex- 
cess in the near infrar e d bac kground. The main conclusion of 
ISalvaterra fc Ferraral l|2006l ) is that the excess near infrared 
background cannot be mainly coming from high-z galaxies 
at z > 7, as they do not see these extremely bright galaxies 
in the NICMOS UDF, ZEN, or Spitzer counts. 

However, their argument does not rule out the 
high-2 galaxies being the origin of the near infrared 
background. Using a simple argument based upon en- 
ergy conservation, we have shown in the previous paper 
(jFernandez fc Komatsull2006D that the near infrared back- 
ground measures only the total light integrated over time, 
and thus one can obtain the same amount of near infrared 
background by having either (i) extremely bright sources 
ov er an extremely short time period, such as those invoked 
bv lSalvaterra fc Ferraral (|2006l ). or (ii) much dimmer sources 
over a much longer time period. While [Salvaterra fc Ferraral 
(2006) have successfully shown that the first possibility is 
ruled out, they have not ruled out the second possibility 
yet. 



2 Q„/(n6r,) ~ 10. 
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5.3 iLe Delliou et al.l (120061 1 

iLe Delliou et all l|2006h predict the luminosity functions of 
Lya emitters at redshifts from 3 < 2: < 6.6, using cosmolog- 
ical simulations coupled with a semi-analytical galaxy for- 
mation model. 

Similar to ours and the other work, they assume that 
the escape fraction of Lya photons are independent of halo 
mass, and find its value, a^sc ~ 0.02, by fitting the ob- 
served luminosity fun ction of Lya emitters at z ~ 3. (See 
iKobavashi et al.]|2007l . for a criticism on this assumption.) 

Since the halo mass function is also an essential ingre- 
dient in the semi-analytical galaxy formation model, and 
they make the same assumption about the escape fraction of 
Lya photons, we expect our predictions and theirs to agree 
well for the same set of parameters. We find that we can fit 
the bright-end of their predicted luminosity functions (their 
Fig. 1 for z — 7) with a population of starburst galaxies, 
Mh/Lboi ~ 1, which is a very reasonable result. 

We believe that our simple model captures the basic 
physics that goes into their model, which is more sophisti- 
cated and complex. 

At a fainter end, however, their luminosity function fiat- 
tens out and our calculations always over-predict the num- 
ber of sources. This is likely due to our assumption about a 
constant mass-to-light ratio. It is expected that this assump- 
tion breaks down once a large mass range is included in the 
analysis. The most economical way to improve our model 
is to introduce a second free parameter, a slope of mass-to- 
light, such that L oc M'^ , for instance. As the observations 
improve in the future, a two-parameter model such as this 
should be used. 



6 CONCLUSIONS 

A simple model based upon the halo mass function coupled 
with a constant mass-to-light ratio fits the luminosity func- 
tions measured and constrained by the current generation 
of narrow-band Lya surveys from 5.7 ^ 2: ^ 8.8. We have 
explored various metallicities and stellar mass spectra. 

The inferred mass-to-light ratios are consistent with no 
evolution in the properties of Lya emitters or opacity in the 
IGM from 5.7 ^ z ^ 7. Therefore, the current data of the 
luminosity functions do not provide evidence for the end of 
reionization. The data at z = 8.8 do not yield a significant 
constraint yet. 

These mass-to- light ratios suggest that the Lya emitters 
discovered in the current surveys are either starburst galax- 
ies with only a smaller fraction of Lya photons escaped from 



galaxies themselves and the IGM, aescf 



1/7 



' 0.01-0.05, or 



normal populations with a fair fraction of Lya photons es- 



caped, aescE 



1/7 



0.5 — 1. The luminosity function alone 



cannot distinguish between these two possibilities. 

For the duty cycle of order unity, e ~ 1, the observed 
equivalent width of Lya line indicates that starburst popu- 
lations are consistent with low metallicity populations with 
Z < 1/50 Zq, while normal populations are consistent with 
high metallicity populations. The other solutions are still al- 
lowed for e < 1. Note that a rece nt study of the SED of Lya 
emitters bv lNilsson et al.l l|2007l ) shows that the Lya emit- 
ters a,t z = 3.15 are consistent with a very low metallicity 
population, Z = 1/200 Zq. 



To constrain the properties of Lya emitters fur- 
ther, one should use asymmetric absorption features of 
the measured Lya line profiles to distinguish between 
them (Miralda-Escuddll998l: [j^alda-fe^ fc Rees, .19981 : 
lHaimanll2002l : ISantosll2004l : lTasitsiomill2006l '). The best way 
to break degeneracy between e and Oesc is to detect the de- 
viation of the cumulative luminosity function from a pure 
power-law. In order to do this it is crucial to determine the 
bright end of luminosity function more accurately. 

We disagree w ith the conclusion reached by 
ISalvaterra fc Ferraral (|2006l ) that no detection of Lya 
emitters at z = 8.8 excludes the excess near infrared 
background being produced by galaxies at z > 7. While 
they have excluded the excess background coming from 
extremely bright starburst galaxies with Mh/Lboi ~ 10~^ 
and the lifetime of 2 Myr, their argument does not exclude 
another possibility that the excess background originates 
from galaxies with Mh/Lboi ~ 0.1 — 1 and the lifetime 
comparable to the age of the universe ai z > 7. As the 
near infrared background measures only the total amount 
of light integrated over time, both scenarios result in the 
same amount of background light. As we have shown in this 
paper, the latter scenario is consistent with all the existing 
Lya surveys from 5.7 ^ « ^ 8.8. 

There is a subtle hint that 20 — 30% more Lya pho- 
tons survived from z — 6.5 than from z = 5.7. A number of 
factors need to be checked carefully before this conclusion 
is taken seriously: the completeness correction and spectro- 
scopic confirmation rate of the observed luminosity function, 
the shape of Lya line profiles (which is however not quite 
enough to make the effect go away; see Appendix|X]| , and ac- 
curacy of the evolution of the theoretical halo mass function 
in these redshifts. In addition, more elaborated theoretical 
models such as those described in Sec. 12.11 may be neces- 
sary to test reality of this effect, while it is interesting that 
the model with a duty cycle has also shown a similar trend 
IStark et al.ll2007l '). 

Our method should provide a simple tool for interpret- 
ing the galaxy number count data in terms of the mass-to- 
light ratio. Or, for a given mass-to-light ratio as constrained 
by the existing data, it can also be used to predict the lu- 
minosity functions, and thus it helps to design future Lya 
surveys such as a follow-up survey at 2: = 7 with Subaru, at 
z — 8.8 with VLT, and a new survey at even hi gher redshifts 
with the James Webb Space Telescope (see also lBarton et all 
|2004 for an alternative way of making forecasts). 
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APPENDIX A: ON THE PROFILE OF LYa 
LINES 

Throughout this pa per we have adopted a theoretical hne 
profile computed bylLoeb fc RvbickH (|l999l ) and later fit by 
ISantos et al.l (|2002l ') (their equation 15). This line profile is 
fairly broad (see dot-dot-dot-dashed lines in Figure U). The 
physical origin of this broadening is a combination of the 
IGM scattering and cosmological redshift. Since their un- 
derlying assumption that the IGM around sources is neutral 
may not be always valid, a care must be taken when one uses 
their theoretical profiles. The broadening is reduced signifi- 
cantly when the IGM around sources is ionized. 

The shape of the line profile affects our analysis through 
the bandwidths of instruments. Since we are dealing with 
narrow-band filters, instruments miss a large fraction of Lya 
photons if a line profile is broader than their bandwidths. 
Therefore, if we assumed erroneously that the line profile 
was too broad, then the inferred M^/Ltoi would be too low. 

In order to quantify a possible uncertainty regarding 
the shape of the line profile, we explore the extreme case 
where a line profile is a delta function at 10.2 eV. Since the 
instruments would not miss any Lya photons, the luminosity 
within the band, Lband, would increase. 

Table lAll shows that a delta-function line profile in- 
creases the luminosity within instrument's bandwidths sub- 
stantially. For metal-free stars with a heavy mass spectrum, 
the luminosity increases by a factor of 2 for Subaru, a factor 
of 3 for LALA, and a factor of 5.5 for ISAAC. 

In Table IA2I we report the inferred 
{Mh/Lboi)iceesc€^^'')^^ from assuming a delta function 
line profile. These values should be compared with those 
in Table [5l While there are changes in the inferred 
{Mh/ Lboi)iciesce^^'')~^ at the level of a factor of a few, the 
main result from our analysis does not change: the Lya 
emitters discovered in these surveys are either starburst 
galaxies with aeac^^^"' ~ 0.01 — 0.1, or normal galaxies with 
a^sce"-^'' ~ 0.5 - 1. 

On the other hand, a hint that (Mh/ Lboi){cteac(-^^~')^^ 
at z = 6.56 is smaller than that at z — 5.7 is now less 
significant: it's only a 10-20% effect, rather than a 20-30% 
effect. 
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Table Al. The same as Table |4] but for a Ly« line profile being a delta function at 10.2 eV. Small differences between 
observations are due to differences in the continuum flux within bandwidths. 
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Salpeter 
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Table A2. The same as Table |5] but for a Lya line profile being a delta function at 10.2 eV. 
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